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Modified f(R) gravity is one of the most promising candidates for dark energy, and even for 
the unification of the whole cosmological evolution, including the inflationary phase. Inside this 
class of theories, the so-called viable modified gravities represent realistic theories that are capable 
to reproduce late-time acceleration, and satisfy strong constraints at local scales, where General 
Relativity is recovered. In the present work, the cosmological evolution for some of these models 
is analyzed, which indicates that these f(R) theories may lead to a phantom phase in the universe 
evolution. Furthermore, the scalar-tensor equivalence of f(R) gravity is also considered, which can 
provide important properties on this kind of models. Moreover, the possibility of the occurrence of 
future singularities as well as the so-called Little Rip are studied. 
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I. INTRODUCTION 



m 

Since the discovery by two independent groups of the accelerating expansion of the universe through observations 
of the luminosity distance of Supernovae IA, as well as by some other independent observations (such as the Cosmic 
\ Microwave anisotropics), a great effort has been made to understand what kind of mechanism is producing this 
CJ^ 1 anomalous behavior on the expansion of the universe, encompassing the mystery under the name of dark energy. 

Since the value of the observed dark energy density is so large in comparison with the one predicted by quantum field 
bJQ" theories for the vacuum energy density, it is quite reasonable to think that an unknown contribution is affecting the 
l— ~~ cosmological expansion, either by a new field such as quintessence models or vector fields, or by some kind of modifica- 
I , tion of General Relativity (GR) which has effects only on cosmological scales (for a review on modified gravities see |l|). 

>: 

Then, as a consequence, along recent years the study of modified gravity theories has become very popular, and 
■ a major point for part of the scientific community, in order to understand the problem and the possibilities to 
reconstruct a gravitational theory capable to reproduce the dark energy epoch, and even the inflationary phase. In 
. this sense, f(R) theories seem to be the best candidate because of their simplicity, generalizing the Hilbert-Einstein 
t""^ 1 action to a more complex function of the Ricci scalar, and the ease of reconstructing the corresponding action that 
reproduces a particular cosmological evolution (see Refs.0-11])- Some of these models can reproduce exactly the 
same behavior as ACDM model (see Ref. Q), or mimic a cosmological constant at the current epoch, unifying also 
the inflationary phase (see 0, 0])- Furthermore, as suggested by observations, the possibility that the dark energy 
behaves as a phantom fluid, i.e. its equation of state (EoS) is less than -1 right now or in the near future, is not 
excluded and has been widely explored (see Ref. [8, 9]). In such a case, the universe evolution may enter in a phantom 
phase, producing a super-accelerating expansion that may lead to some kind of future singularity, as the so-called 
Big Rip (for a classification of future singularities, see Ref. [H|). In this sense, f(R) models that cross the phantom 
barrier and lead to future singularities have been also studied in Ref. 11]. 



Nevertheless, the high precision of local gravity tests makes that any modification of GR should pass very 
strong constraints, and recover GR in a certain limit (see Ref. [Hj]). In this sense, some viable f(R) models have 
been proposed, which are capable to reproduce a realistic cosmological evolution, and prevent violations of local 
gravitational tests and instabilities in the presence of matter distributions (see Refs. [13l - [l5| ). Furthermore, some of 
these models are capable to reproduce both the dark energy epoch as well as the inflationary phase (see Ref. 16]). 
However, viable f(R) gravity may lead to a phantom phase in the present/future time, which may produce some 
kind of future singularity (see Ref.[17|). Moreover, in the era of precise cosmology, a deeper study on the large 
number of dark energy models is completely mandatory in order to solve a problem of degeneracy, since most of 
dark energy candidates can reproduce the same cosmological evolution with no significant differences. In the case of 
modified gravity, some studies about the stability of the cosmological solutions [ll| , or the exploration of cosmological 
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perturbations (19| , can provide some restrictions on the form of the gravitational action. 

In the present work, we focus on the study of the so-called viable f(R) models, exploring their cosmological evolu- 
tion, and its scalar-tensor representation. It is well known that /(i?) gravity can be redefined as a kind of Brans-Dickc 
theory in terms of a scalar field with a null kinetic term (see Refs. [g, 12CJ and references therein), so the phase space 
for Friedmann-Lemaitre-Robertson- Walker (FLRW) metrics can provide some useful information about the stability 
of the cosmological evolution as well as the presence of future singularities in viable f(R) gravity. In this sense, we 
study two characteristic models of the class of viable f(R) theories (see Refs. [Lj|l4(), where the cross of the phantom 
barrier appears as a common factor for both models, and where the possibility of a future singularity is studied. In 
addition, the possibility of the so-called Little Rip is also explored, this is a postulated phase of the universe evolution, 
free of sin gul arities, where a very strong accelerating expansion may lead to break some bounded systems, as planetary 
systems, [21]. Such possibility may introduce some important predictions (and constraints) in the gravitational 
models as it may produce instabilities on theplanets orbits in the near future. The Little Rip has been previously 
considered for several dark energy models [22| . as well as in some contexts of modified gravities, [23| . Here we 
study the possibility of a Little Rip in viable f(R) theories, in particular we study the effects on the Earth-Sun system. 

The paper is organized as follows: section [TT] is devoted to the introduction of viable f(R) theories and their 
main cosmological properties. In section ITOI we explore the evolution of the effective EoS parameter of the universe, 
showing that in general the phantom barrier is crossed. Section llVI deals with the analogous scalar-tensor theory for 
f(R) gravity and the behavior of the phase space. Finally, in section |Vj we discuss the possibility of the presence 
of cosmological singularities according to the results found in previous sections, as well as the occurrence of a Little 
Rip along the cosmological evolution, where we study the effects of the expansion on the gravitational coupling in the 
Earth-Sun system. 



II. VIABLE f(R) THEORIES 

Let's start by reviewing reliable cosmological models for the kind of modified gravities called viable f(R) theories. 
We consider the action corresponding to one of these theories which, aside from the gravity part, also contains a 
matter contribution, namely 



S = J d^x^g [f(R) + 2n 2 C m ] 



(1) 



Here the coupling constant is given by k 2 = 8irG, while L m stands for the Lagrangian corresponding to the matter 
that fills the particular system to be studied. Note that the Hilbert-Einstein action is recovered for f(R) = R in (JTJ|, 
The field equations corresponding to action ([l} are obtained by the variation of this action with respect to the metric 
tensor g^ v , what yields 



R^f R (R) - \g^f{R) 



Uf R {R) - V M V„/*(fl) = n 2 T^ 



Here the subscript r denotes derivatives with respect to R. We assume a flat FLRW metric, 

3 

ds 2 = -dt 2 + a(t) 2 V dx ir 



i=l 



(2) 



(3) 



Then, for this metric, modified FLRW equations are obtained through the 00 and ij components of the field equations 



1 

3?r 



K Pm 



RJr - 1 



SHRJrr 



-3H 2 - 2H = — 

JR 



K 2 p„ 



R 2 f R RR + 2HRf RR 



RfRR + ^(f-RfR) 



(4) 



where dots denote derivatives with respect the time, Hubble parameter is defined as usual H(t) = a/a, and Ricci 
scalar for the metric ([3]) is R — 6 (2H 2 + H). Here, we are interested in studying a subclass of modified gravities, the 
so-called viable f(R) gravities, as they accomplish some indispensable conditions to be considered realistic candidates 
to describe the universe evolution. This kind of f(R) gravities are usually described by actions of the type, 



f(R) = R + F(R) 



(5) 
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This action basically represents Hilbert-Einstein action plus an additional term that should have only effects at 
cosmological scales, while at local scales, General Relativity should be recovered. In order to avoid serious problems 
with known physics, one has to choose the F(R) function such that the theory contains flat solutions and passes also 
local gravity tests (see [IH). Viable gravities are able to satisfy these constraints, and to avoid large instabilities in 
the presence of matter distributions, as well as the anti-gravity regime by imposi ng p ositivity in the first derivative of 
the action, 1 + F R > 0, and the appropriate form of the function F(R) (see [13(-[ia|). Moreover, on the cosmological 
level, these theories are capable to reproduce the effects of dark energy, and even inflation. Note that both equations 
in are written in such a way that -F(i?)-terms are put on the matter side (see Ref. [6|); thus we may define an 
energy density for the extra F(i?)-terms. Hence, the first FLRW equation in Q can be rewritten for the kind of 
actions ([5]) in terms of the cosmological parameters f2j, 

1 = Q m +CIf(r) j where 



n m = , n F(R) = — — | MRFrr - 3H 2 F R . (6) 



1 ( RF R — F . >f;r>r , 

K V 

Then, the first Friedmann equation Q takes a simple form, with two fluids contributing to the scale factor dynamics. 
In addition, the continuity equation V M T Ml/ = for a perfect fluid with an EoS p m = w m p m yields, 

p m + 3H(1 + W m )Prn = . (7) 

Hence, for a particular F(R), the corresponding cosmological evolution can be obtained through equations (|4|), and 
the continuity equation (|7|). 

In absence of any kind of matter, dynamics of the universe are carried out by the i 7 '(7?)-component, which may be 
chosen so that reproduces (or at least contributes) to the early inflationary phase and the late-time accelerating epoch. 
To reproduce the whole history of the universe, the following conditions on the F(R) function have been proposed 
(see P3): 

i ) Inflation occurs under one of the following conditions: 

lim F(R) -» -Aj or lim F(R) -> aR n . (8) 

R— >oo R— >oo 

In the first situation in (J8J, F(R) function behaves as an effective cosmological constant at early times, while 
the second condition yields also an accelerating expansion where the scale factor behaves as a(t) ~ t m , with 



m = 2n'+l-3n 
2—n 



(see Ref. Q). 



ii ) In order to reproduce late-time acceleration, we can impose a similar condition on F(R). In this case, Ricci 
scalar has a finite value, which is assumed to be the current one, so that the condition is expressed as 

F(R ) = -2R F'(R )~Q. (9) 

This basically means that extra terms in the action behave as an effective cosmological constant at the present 
time, p ~ —p. However, as the effective fluid, coming from extra terms in the action, is not exactly a cosmological 
constant, the evolution would be quite different, since the universe may enter in a phantom phase (l7| . or may 
produce oscillations along the cosmological evolution [24[ . We will explore the evolution, and the crossing of the 
phantom barrier in viable f(R) gravity in next sections. 

Hence, under these circumstances, the F(R) term is able to reproduce both accelerating epochs of the universe 
evolution. Here we will focus on the study of two models of this kind of viable gravity, proposed by Hu and Sawicki 
in Ref. [l3j], and Nojiri and Odintsov in Ref. [l4j, whose actions are given by, 

C2{R/Rhs) + 1 l + cR n 

Here, {ci, C2, n} are free parameters and Rhs = ^ 2 p m according to Ref. |13|. while {a, b, c, n} are free parameters for 
the second model in ([TU]) . The first model in (fTU)) has been studied in Ref. Q, where it was proven that the universe 
evolution goes through two different De Sitter points, being one of them stable and the other one unstable, which can 
be identified to the current accelerated era and to the inflationary epoch, respectively. Similarly, the second model 
in (|10p was studied in Ref. [16], also with the presence of an extra field, and was shown that can well reproduce 
both accelerated epochs, mimicking a cosmological constant at the present time. Thus, the models (fTU|) are capable 
to reproduce both accelerated epochs, but a different cosmological evolution than ACDM model is reproduced. In 
next sections, we will study the evolution of both models, specially during the dark energy epoch, focusing on the 
possibility of the crossing of the phantom barrier. We will also study the behavior of their scalar-tensor equivalent, 
where the phase space will be explored. 
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III. COSMOLOGICAL EVOLUTION IN VIABLE f(R) GRAVITY 

In this section, we explore the cosmological evolution for the models considered above. For that reason, it is 
convenient to perform a change of variables in the equations Q, considering the redshift z as the independent 
variable instead of the cosmological time t, where the redshift is defined as usual, 

l + z = ^- 7 (11) 



and ao is the value of the scale factor at the present time t , such that the current epoch corresponds to z — Then, 



the time derivative is transformed as 4l = —(l + z)H-4-, and the first FLRW equation in (j4|) and the continuity 



equation ([7]) yield. 

U 2 p m {z) + R{Z) { R - f + 3(1 + z)H 2 f RR R\z) 



(12) 



(1 + z)p' m {z) - 3(1 + w m )p m (z) = , (13) 

where now primes denote derivatives with respect to the redshift. Then, the Ricci scalar can be rewritten as R = 
6 [2H 2 (z) - (1 + z)H(z)H'(z)] , while the equation (T3]) can be easily solved for a constant EoS parameter w m , 

p(z) = po(l + z)^ 1+w ^ . (14) 

Here po is the value of the matter energy density at the present epoch z = 0, which can be rewritten as op = AsHq f2^. 
Then, we can fit the current values of the cosmological parameters using the observational data 26] , where Hq = 
100 h km s _1 Mpc" 1 with h = 0.71 ± 0.03 and the matter density, fl^ = 0.27 ± 0.04, while the matter fluid is 
considered pressureless (cold dark matter and baryons), such that w m = 0. The equation (fT2")l is a second order 
differential equation on H(z), so by fixing the initial conditions, the corresponding cosmological evolution can be 
obtained through the Hubble parameter in terms of the redshift for both models considered in previous section, and 
explore how the future evolution (with — 1 < z < 0) might be. Here, we explore the evolution of the effective equation 
of state parameter w e f / , which is defined as 

_ PF(R)+Pm _ 2H(t) _ 2(1 + z)H'(z) 

eff ~ PF lR )+Pm~ 3H*(t)- 1+ 3H(z) ■ (15) 

For simplicity we redefine the Hubble parameter as H(z) = H h(z), such that h(0) = 1, while the initial condition 
on the first derivative h'(0) can be fixed by extrapolating a particular model at z = 0. In this sense, if we assume 
that ACDM is a good description at z = 0, FLRW equations are given by, 

k 2 A k 2 

H 2 = —Pm + 3 , H= - — p m ■ (16) 

Performing the change of variables t — > z, and H(z) = H h(z), the second equation in ([T^| at z = yields, 



We are interested to explore the cosmological evolution for the models (flOf by assuming that the universe mimics 
ACDM model (fTTH) at present time (z — 0), so that the initial conditions (jTT)) are assumed. Moreover, we may 
compare such assumption with the one that supposes the universe has already crossed (or will do very soon) the 
phantom barrier, i.e. the universe is already dominated by a kind of phantom fluid at z = 0, a possibility allowed 
by observations. In such a case, the first derivative of the Hubble parameter with respect to the redshift should be 
approximately zero or negative, h!(fi) < 0. 

Let's now analyze the Fhs{R) model in (|10[) . assuming c\ — 2 and C2 = 1, which are dimensionless parameters, and 
a power of n = 1. In Fig. 1, the evolution of the EoS parameter (fT"5)) is plotted with respect to the redshift for initial 
conditions that mimic ACDM model at z = 0. In spite of the evolution of the EoS parameter is quite similar to the 
ACDM model for positive redshifts, the universe enters in a phantom phase in the future (negative redshifts), where 
the EoS parameter turns out less than — 1, and it oscillates close to z — — 1 as shown in fig. Ilbl in more detail. In 
Fig. 2, the EoS parameter (| 1 5|) is also evaluated, but assuming different initial conditions at z = 0. In this case, the 



(a) (b) 

Figure 1: Evolution of the EoS parameter (1 1 5 1> as a function of the redshift z for the model Fhs{R) where n = 1, ci = 2 , cz = 1, and 
initial conditions h(0) = 1 and h'(0) = 0.4, according to 1171 . Figure Hal shows the evolution from the past to the future, where the 
phantom barrier is crossed for negative redshifts. In fig. Ilbl we can see in more detail the range —1 < z < —0.8, where the phantom 

barrier is clearly crossed, and the EoS oscillates. 




Figure 2: Evolution of if e // for the model Fhs(R), where n = 1, ci = 2 , ci = 1, and assuming that the universe has already entered 
in the phantom phase, h(0) = 1 and h'(0) = —0.1. As shown in figure [2al the EoS crosses the phantom barrier even before 2 = 0, and it 
oscillates along the barrier in the future. In more detail, we can see the future evolution in fig. I2bl which ends outside of the phantom 

region 



universe has already entered in the phantom phase at z = 0. The evolution for the second case, plotted in Fig. 2, 
shows that the amplitude of the oscillations start before, and they are larger than in Fig. 1, while the evolution for 
redshifts z > 0.5 is quite similar. Hence, both examples can give accurate descriptions of the universe evolution in 
the past, but they predict clearly different behaviors for the future. Furthermore, both examples present oscillations 
around the phantom barrier when the redshift is close to z = — 1, as shown by Figs. llbH2bl 

In Figs. 4 and 5, the evolution of the EoS parameter for the NO model (I10|) is evaluated. As in HS model, firstly 
we assume initial conditions given by h'(0) — (Fig. 3), and then, in Fig. 4, h'(0) ~ —0.1. For NO model, not all 
the free parameters are dimensionless as above, but {a, c} have dimensions of H~ 2 while b is effectively dimensionless. 
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(a) (b) 

Figure 3: Evolution of the EoS parameter for Fmo(R)- We have assumed n = 1, a = 0.1/Hq , b = 1, c = 0.05/-Hq and initial conditions 
similar to ACDM model at z = 0, ^(0) = 1 and h'(0) = 0.4. As shown, this model behaves very similarly to Hu-Sawicki model, fig. ?? is 
similar to fig. Hal while the graphic l3bl shows that the future evolution enters in the phantom region, and bounces above the barrier, with 

no oscillations. 




As shown in Fig. 3, the behavior of NO model is similar to the HS viable model: the EoS evolutes similar as ACDM 
model for z > 0, while it enters in the phantom region when z < where main differences with respect the above 
model are found. According to the Fig. I3b| the phantom barrier is crossed, but then, it bounces above w = — 1, and 
no oscillations are produced as in Fig. Ilbl For the second case (Figs. 4), where the crossing of the phantom barrier 
is assumed to occur before z = 0, the EoS parameter for NO model oscillates, and the universe crosses the phantom 
barrier in positive redshifts, ending above the barrier w = — 1, similarly as HS model, shown in Fig. I2al 

Figures 5 and 6 show the evolution of the cosmological parameters {f2 m ,f2j?} with respect to the redshift for the 
Hu-Sawicki and Nojiri-Odintsov models respectively for the different initial conditions considered above. Note that 
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Figure 5: Evolution of cosmological parameters {Q m , £2f} defined in (O as a function of the redshift z for the HS model with n = 1, 
c\ = 2 , C2 = 1, and assuming initial conditions, 7l(0) = 1 and h'(0) = 0.4 (fig. |5al , and h(0) = 1 and h'(0) = —0.1 (fig5b). Both 
evolutions seem very similar, but in the right one, the cosmological parameter Qp decays for < z < 1, but it starts to increase again for 

z > 1, clearly in contradiction with standard model. 





-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 
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(b) 



Fi gure 6i Evolution of cosmological parameters {r2 m (z), SJf (2)} for Fno(R) with n = 1, a = 0.1/Hq , b = 1, c = 0.05///q, and 
assuming initial conditions, /i(0) = 1 and fo'(O) = 0.4 (fig. I6a[ l. and /i(0) = 1 and h'(0) = —0.1 (fig. I6b[l . In fig. I6bl the transition to the 
dark energy epoch occurs before than in fig. I6al a consequence of the initial conditions. 



the choice of the initial conditions affect the cosmological evolution for both models, while for the HS model, both 
cases present the same evolution for negative redshifts, for positive ones, the evolution is different as shown in Fig. l5b| 
where the cosmological parameters seem to oscillate when, z > 0. In the NO model, when phantom initial conditions 
are assumed, tip tends to dominate completely the universe, while for positive redshifts, it becomes negative, fig. I6bl 

Note that, as shown in figures 1-4, viable modified gravities ([10]) produce some oscillations along the cosmological 
evolution, a fact pointed out before for this kind of models in Ref. [24| • In general, we have shown that the transition 
to the phantom epoch occurs. Moreover, the election of phantom conditions at z = gives an anomalous behavior of 
the cosmological parameters for both models when z > 0. 
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IV. SCALAR-TENSOR REPRESENTATION OF f(R) GRAVITY 

It is well known that f(R) gravity is equivalent to a kind of Brans-Dicke theory with a null kinetic term, and a 
scalar potential, (see for example, [y,[2(| and references therein), 

^V^g [<t> R ~ V(4>) + 2n 2 C m ] . (18) 



By varying the action with respect to the scalar field <fi, the relation between both theories is obtained, and f(R) 
action |T]) is recovered, 

R = V'(4>) -> <j) = <j){R) , => f(R) = <j){R)R - V{<t>{R)) . (19) 

While the scalar field and the potential are related with the particular f(R) action by, 

= f R (R) , V(<f>(R)) = fn(R) R - f(R) . (20) 

Then, by a particular f(R) action, we can reconstruct its equivalent scalar-tensor theory Q18]l . so that we can study 
the viable models studied above by analyzing their scalar-tensor counterpart. Moreover, by varying the action (|18j) 
with respect to the metric g^ Vl field equations are obtained, 



R 



3D4> = K 2 r< m ' + - 2V{<(>) , (21) 

dep 

where the second equation is the trace of the field equations in (|2"Tj) . Hence, for a fiat FLRW metric ©, the Friedmann 
equations yield, 

3H 2 = ^-( K 2 p m -3Hj ) +^V(4>)) , -3H 2 - 2H = I L 2 Pm + '<j> + 2H0 - \v{<t>)\ . (22) 



while the trace equation in (|21l) is given by, 

30 = K 2 ( Pm - Pm ) - <W' + 2V~ 9H<j> . (23) 

Let's reconstruct the corresponding scalar-tensor theory for the particular models considered in the previous section. 
By the expressions (f2"0"l) . the corresponding relation <p(R) for the models (fit)]) are given by, 

C1C2R C\ cR(aR—b) aR aR—b 

where we have assumed a power of n = 1 for both models (fTQ)) . Then, the scalar potentials ([20]) yield, 



feW = i±^M^MEl %s , ^ oW ^ 2a + c(1 + b ^ )±2 2 ^ a+6c ) (25) 

Hence, the scalar-tensor representation for both models is not uniquely defined, but the potentials exhibit two branches, 
which in principle do not affect the cosmological evolution, but it may influence the behavior of the phase space. In 
addition, both models introduce a boundary condition on the value of the scalar field, being 4> < 1 for the HS model, 
and 4> < a + c / c f° r the NO model, a limit where both branches of the potentials converge, as can be seen in Fig. 7. 

Moreover, the evolution of the scalar field as a function of the redshift, <j>(z) — fii(R(z)), is the same regardless of 
the sign in the potentials (|2"5|) . as shown in Fig. 8. Note that there is a direct correlation between the behavior of 
the scalar field and the evolution of the EoS parameter as can be shown comparing figs. 1 and 3 with fig. 8. For the 
HS model, the EoS parameter (Fig. l)presents some oscillations, specially for negative redshifts around the phantom 
barrier, while the evolution of the scalar field also oscillates according to Fig. [5a] In the NO model, the linear behavior 
of the EoS parameter observed in Fig. 3 is also translated in a linear evolution of the scalar field with respect to the 
redshift, see fig. I8bl Let us now analyze the phase space {H, <fi} for this kind of viable models in vacuum, where 
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Fi gure 7'. Scalar potentials V((f>) in terms of Hq, of the HS model with n = 1, c\ = 2 , C2 — 1 (figure I7a|l . and of the NO model with 
n = 1, a = 0.1/ HS , b = 1, c = 0.05/ HS, fig, I7bl Both potentials are not uniquely defined, where each branch behaves very different, 

converging to the boundary of the scalar field <f>. 




(a) (b) 

Figure 8: Evolution of the scalar field (p(z) for the HS model with n = 1, ci = 2 , c% = 1 (a), and for the NO model (b) with n = 1, 
a = 0.1/Hq , 6 = 1, c = 0.05/Hq. In both cases, the initial conditions are: /i(0) = 1, and h'(0) = 3f2^/2. For each example, we can 
easily see the correlation between the evolution of the EoS parameter in Figs. 1 and 3, and the evolution of the scalar field for each 

model respectively 



both branches of the scalar potentials are considered. For simplicity, we assume vacuum FLRW equations, so that by 
combining equations (f22]) - (f23|) , we get, 



H = -2H 2 + \V{<t>) , 4> = ^~ 



G 



3H 



-3H 2 0+lv(<P) 



(26) 



Hence, by combining both equations in (|26|) . 



dH 



6 



2H 2 



(27) 



3H 
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(a) (b) 



Figure 9: Phase space H((f>) for the HS model, where different initial conditions are considered, and both branches of the potential are 
studied. The left figure corresponds to V+(<fi) and the right one to V— (cj>). In both cases, we have considered initial conditions close to 
the boundary of the scalar field: (a) #(0.95) = #g, (6) #(0.95) = WH$, (c) #(0.95) = 100#g. 




(a) (b) 

Figure 10: Phase space H(4>) for the NO model. As above, different initial conditions are considered, and both branches of the 
potential are studied. The left figure corresponds to V+(<j>) and the right one to V— (</>). We have considered initial conditions close to the 
boundary of the scalar field: (a) #(2.95) = H$, (6) #(2.95) = 10#g, (c) #(2.95) = I00#^. 

This equation describes the phase space for a particular scalar potential, and provide the information about the 
behavior and stability of a particular model. Specifically, we can see in Fig. 9 that regardless of the choice of initial 
conditions, the behavior of the Hubble parameter tends to a stable point, close to zero. However, as shown above, the 
scalar field takes values very distant of the stable point, so the sensibility to the initial conditions can be important 
along the cosmological evolution naturally, although the evolution presents the same tendency independently of 
the initial conditions. In fact, we can observe some small fluctuations along the evolution of H with respect to 
<j>. Moreoever, both branches of the potential Vhs in (E3 may present different convergence points in spite of the 
differences may not be detected in the range of values that the scalar field takes along the cosmological history(see 
Fig. 8). For the NO model in Fig. 10, a similar behavior is found, where different convergence points are found for 
V±, as well as larger fluctuations for V+, while the behavior is similar in both cases, and independently on the initial 
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conditions (assumed close to the boundary of the scalar field). 

Hence, the cosmological evolution in viable f(R) gravity shows a similar behavior regardless of the model, and even 
of the initial conditions, where the crossing of the phantom barrier is a common element in all of them, but with 
no sign of possible future singularities, such that viable f(R) gravity represents a serious candidate for explaining 
the dark energy epoch. The tiny differences between models seem very clear when the future evolution of the EoS 
is studied, and whose phase space does not seem to provide useful information for its discrimination. However, the 
stability shown in the phase space for both models suggest that no future singularity will occur in the future. In next 
section, we discuss this point in more detail, and we analyze the possibility of the occurrence of a Little Rip. 



V. FUTURE SINGULARITIES AND LITTLE RIP IN VIABLE f{R) GRAVITY 



In recent years, the study of future singularities has become a major task since a lot of cosmological models, capable 
to describe a realistic cosmological evolution satisfying observational constraints, lead to the so-called phantom phase, 
w < — 1, and may end in some type of future singularity (see Refs. Q). A classification of future singularities was 
presented in Ref. [lOj . 

• Type I ("Big Rip"): For t — > t s , a — > oo and p — > oo, \p\ — > oo. 

• Type II ("Sudden"): For t — > t s , a — > a s and p — > p s , \p\ — > oo. 

• Type III: For t — > t s , a — > a s and p — » oo, \p\ —> oo. 

• Type IV: For t — > t s , a — > a s and p — > p s , p — » p s but higher derivatives of Hubble parameter diverge. 

Future singularities have been widely studied, in particular, they may occur for some f(R) models (see Refs. @,|Tl[). 
Nevertheless, the fact that the universe may cross the phantom barrier at w = — 1 is not sufficient condition to present 
a future singularity (see Ref. Q). In section Hill we have shown that in general, viable models of f(R) gravity leads 
to a phantom epoch in the future, where the EoS parameter becomes less than -1, but no future singularity occurs 
since Hubble parameter and its first derivative remain finite, or in other words, according to the classification, there 
will not be a Big Rip, Sudden singularity or Type III singularity for the viable models studied above. In addition, 
type IV singularity is not completely excluded, but this kind of singularity, that affects to higher derivatives of the 
Hubble parameter, may not have any physical consequence, at least observationally as pointed out in Ref. (Tfj, so its 
importance in comparison with the others is lower. 



Nevertheless, the absence of future singularities for a particular model is not enough to avoid some critical conse- 



quences as the so-called Little Rip. The Little Rip is a postulated phase of the universe evolution Ref. [21], when a 
very strong accelerating expansion would lead to break some bounded systems, as the Solar System or the galaxies, 
but where all physical magnitudes remain finite (see also Ref. [22]). This special scenario has been also explored in 
modified gravity, where the occurrence of this scenario seems very possible for some particular models (see Refs. 23]). 
Here we are interested to study the possibility of a Little Rip in viable f(R) theories, in particular in the models 
studied in previous sections. Firstly, we should characterized in some way the strength of the force of the cosmological 
expansion in order to compare with those forces of a bounded system, as for example the newtonian force in Solar 
System. In this sense, we may write the force that a body of mass m may feel due to the cosmological expansion by 
the relative acceleration between two points separated a comoving distance r, 



(H 2 +ii} . (28) 



Then, we may compare this force with the bounding force of a particular system. For example, in the case of the 
Solar system, the force that maintains the coupling between Sun and a body of mass m can be characterized by the 
Newtonian law, 

F * = G ^P- ■ ( 29 ) 

where Mq = 1.9 10 30 kg is the mass of the Sun, and r is the distance between the Sun and the body m. Here, in both 

cases it is assumed that Fi = \Fi\, as both forces have radial directions, but opposite sign. Hence, we may compare 
both quantities along the cosmological evolution, 

^rf = c^Hl [h(zf - (1 + z)h(z) h'(z)] , (30) 
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where recall that we defined H(z) — Hq h(z) in section Hill where Hq = 100 h km s -1 Mpc -1 with h = 0.71 ± 0.03. 
Let's compare both forces for the Earth-Sun system, where r = 149.6 x 10 9 m is the Earth-Sun distance, at the 
present time, 2 = 0, 

%^ = 1.4 x HT 22 [MO) 2 - h(0)ti(0)] , (31) 
Then, by assuming the initial conditions from section ITTll h(0) = 1 and a) h'(0) = §0^ , and 6) ft'(0) = —0.1, 

(i) %HL ~ 0.82 x 10- 22 , (it) %^ ~ 1.95 x 10" 22 . (32) 
Hence, both values are too small to produce any kind of instability on the orbits of the planets, and in particular on 




(a) (b) 



Figure 11: Evolution of the ratio F "p^ m for the Earth-Sun system in the HS model, where h(0) = 1 and: (a) h'(0) = |f2m , (b) 

h'(0) = -0.1. 

the Earth orbit. In addition, the future cosmological evolution does not seem to affect significantly to the Earth-Sun 
system whatever the f(R) model is used, and independently on the initial conditions imposed at z — as shown in 
Fig. 11 (Hu-Sawicki model) and Fig. 12 (Nojiri-Odintsov model). Nevertheless, some tiny differences are observed 
between both models. For the HS model, the ratio ([30| behaves as a damped oscillator, being the amplitude reduced 
in the future (negative redshifts) , a non expected behavior according to the evolution of the EoS parameter shown in 
section Hill However, the amplitude is too small to affect the Newtonian force between the Earth and Sun. Moreover, 
assuming that the universe has already entered the phantom phase at z = affects directly to the evolution of the 
ratio (I30[) . where the force of the expansion is larger near z = 0, and even negative in the past, where it is assumed that 
the dust matter dominates. In the case of the NO model, different initial conditions affect directly to the evolution 
of the expansion force, being maximum in both cases near 2 = and decreasing when the universe approximates 
to z = — 1. In both cases of the NO model, the force becomes negative around 2 = 1, where dust matter domi- 
nates, but oscillates for small redshifts, where the amplitude is larger in the case of initial phantom conditions, fig. U2bl 

Hence, viable f{R) gravity does not seem to affect the bounding systems, as Solar System. Since the expansion 
force is proportional to the distance between two points, it may affect to some objects far away from the center of 
galaxies as are weakly coupled to them. However, this tiny effect does not seem to be really detectable in comparison 
to other effects. 



VI. CONCLUSIONS 



Along the present paper, we have studied the so-called viable f(R) models by analyzing two characteristic 
examples. We found that the cosmological evolution of the effective EoS for the universe enters a phantom stage 
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in the present or near future, independently of the initial conditions imposed. In this analysis, we have seen that 
both models can reproduce quite well the cosmological evolution, similarly to ACDM model at z = 0, where we have 
imposed some particular conditions in order to fit the model with the observational constraints. As shown in section 
IIII1 when the phantom transition is assumed to occur before z — 0, the behavior of the cosmological parameters 
{il m ,ilp} when z > present some inconsistencies for both models, such that those conditions may be excluded. 
On the other hand, the transition to the phantom phase does not imply directly the occurrence of a future singularity, 
since the Hubble parameter and its first derivatives remain finite. The possibility of a IV type singularity is 
not completely excluded as it implies an infinite number of higher derivatives of the Hubble parameter, but its 
possible physical effects leave unaffected the cosmological evolution. Even in the case that the transition to the 
phantom era has already occurred at z — 0, future evolution seems to be free of singularities. However, the choice 
of the initial conditions at z = may affect the past evolution, as shown in the right fig. I5bl for the HS model, 
where the evolution of the cosmological parameters seems to deviate from the correct evolution for positive redshifts. 

Furthermore, by using the equivalent picture of f(R) gravity as a scalar-tensor theory, we have analyzed the 
evolution of the scalar field in Sect. II VI we note a correlation between the cosmological evolution shown in section Hill 
and the evolution of the scalar field. In addition, we have seen that in general, this kind of viable modified gravities 
are degenerated in their scalar-tensor representation, giving two potentials that impose a bound on the scalar field 
and tends to the same value. Nevertheless, by exploring the phase space in vacuum, we found that independently on 
the initial conditions, both branches tend to stable points, encouraging the idea that no future singularity occurs. 

In the last section, we discussed the possibility of the occurrence of a future singularity, but also the so-called Little 
Rip, an stage of the universe evolution that may affect some bounded systems due to the strength of the expansion, 
but where all the physical magnitudes remain regular. We have explore this possibility by comparing the strength 
felt by a body with respect to another point of the universe, and the newtonian force of a typical planet system. 
For the Earth-Sun system, the expansion force is too small to produce any effect on the orbit. However, as pointed 
above, some systems weakly coupled may have some tiny effects, although the possibility to detect and distinguish 
those effects from others seems not possible. 

Hence, at the present work, we have performed a deep analysis on viable f{R) gravities, showing that the models 
considered here are capable to reproduce the late-time acceleration, entering in a phantom phase in the present or near 
future, a fact that may be a way to distinguish from other dark energy models. Furthermore, the absence of future 
singularities, and instabilities in the orbits of the planets, provide a strong support to the so-called viable modified 
gravity as a serious candidate for explaining the dark energy paradigm. 
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